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Abstract 



We calculate the decay rate of scalar and pseudoscalar Higgs bosons into a pair of 
gluinos, within the Minimal Supersymmetric Standard Model. In the theoretically and 
experimentally allowed light gluino window, m g ~ 3-5 GeV, gluino pairs can completely 
-C dominate the decays of the light scalar Higgs boson and play a prominent role in the 

decay of the pseudoscalar Higgs boson. This would alter the limits obtained from Z 
decays on the lightest CP-even and CP-odd Higgs bosons, and could jeopardize the 
search for these Higgs particles at future hadron colliders. In contrast, the branching 
ratio for the two-body decay of Z bosons into pairs of light gluinos is less than 0.1%. 
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1. Introduction 

The last couple of years have seen renewed interest [TJ in the possible existence of very light 
gluinos, with mass m g < 5 GeV. A recent analysis by Farrar concludes that only the 
upper end of this range is still open, roughly 3 GeV < m g < 5 GeV. However, even this 
window would be sufficient to allow for a substantial modification of the running of as, 
which was the original motivation Q for this scenario. Moreover, since the decays of light 
gluinos might only produce a small amount of missing energy ||, bounds on squark masses 
derived from searches Q for events with large missing Et at pp colliders are not valid in this 
scenario. 

In this paper we point out that the existence of light gluinos could also substantially 
complicate the search for the Higgs bosons predicted by supersymmetric models. To be 
specific, we work within the framework of the minimal supersymmetric standard model or 
MSSM, which contains two SU(2) doublets of Higgs superfields [|], ||. We will not assume 
unification of gaugino masses, which would impose severe constraints on models with light 
gluinos J7|, nor do we require radiative breaking of the SU(2) symmetry, which appears to 
be inconsistent with the existence of light gluinos, at least in the minimal model ||. We find 
that the gg final state can completely dominate the decay of the light scalar Higgs boson; in 
some cases the branching ratios for the usually dominant bb and t + t~ final states are reduced 
by a factor of 1,000 or more. Light gluinos might also be produced in 70% of all decays of the 
pseudoscalar Higgs boson of the MSSM. This has obvious consequences for Higgs searches 
that rely on b or r tagging. The decays of supersymmetric Higgs bosons to gluino pairs have 
been discussed in the past in ref. ||. However, the possibility that the gluinos are light and 
the fact that the mixing between the two supersymmetric scalar partners of the top quark 
(which, as we will see later, is the most important feature in this decay) might be large have 
not been considered. Therefore the Higgs — > gg branching ratios obtained in ref. || were 
very small. 

It has been known for some time |TIJ that light gluinos can also be produced in two- 
body decays of Z bosons; this decay has recently been studied in ref. [fTTJI . We re-compute 
the corresponding partial width, including the effect of mixing between SU(2) doublet and 



singlet squarks. Our numerical results agree with those of ref. |TlJ , but our final expression 
is much more compact. We find a maximal Br(Z — > gg) ~ 6 ■ 10~ 4 within our assumptions; 
such a small branching ratio can only be probed by a dedicated search for gg final states. 
Unfortunately the failure to detect Z — > gg decays would not exclude the possibility that 
gluino pairs contribute significantly or dominantly to Higgs boson decays. 

The remainder of this article is organized as follows. In sec. 2 expressions for the partial 
widths for the decays of Higgs and Z bosons into pairs of light (effectively massless) gluinos 
are given. In sec. 3 we present numerical results for the separate branching ratios and discuss 
possible correlations between them. Sec. 4 contains a brief summary and conclusions. The 
Appendix collects expressions for the decays of Higgs bosons into massive gluinos, as well as 
the necessary two- and three-point functions. 



2. Formalism 

The diagrams contributing to the one-loop induced decays of Higgs and Z bosons into a pair 
of gluinos are shown in Fig. 1. Since the gluinos are of Majorana nature and hence identical 
fermions, one has to antisymmetrize the decay amplitude; this is achieved by adding the 
contributions of the two diagrams where the momenta of the gluinos are interchanged and 
by multiplying these contributions by an overall minus sign. In the case of Higgs bosons, 
this reduces to just multiplying the amplitudes of the diagrams shown in Fig. 1 by a factor 
of two; for the Z boson one needs in addition to discard the terms that are not proportional 
to 75. 

In Fig. 1 one has to include the contributions of both squarks of a given flavor. As well 
known |I2[ , the supersymmetric partners of left- and right-handed massive quarks mix. The 
mass eigenstates q\ and g 2 are related to the current eigenstates qi and q R by 

Qi = h cos 9 q + q R sin 9 q , q 2 = -q L sin 9 q + q R cos 9 q . (1) 

The mixing angle 9 q as well as the masses m qi , m q2 of the physical squarks can be 
calculated from the following mass matrices^, written in the convention of ref . [|l3j : 



kaI = ( m i L + m t + 0-35/^ -m t (A t + fi cot/3) \ 

1 \ -m t (A t + n cot/3) m 2 R +m 2 t +0.16D z )' [ ] 

ka2 _ ( m i L + m b~ 0A2D Z -m b (A b + /itan/3) \ 
M b ~ \ -m b {A b + /itan/3) m\ R + m\ - 0.08£> z J ' ^ ' 

where D z = M z cos2/3, tan (3 being the ratio of the vacuum expectation values of the two 
neutral Higgs fields of the MSSM 0. mj i b are soft breaking masses, A b j are parame- 
ters describing the strength of nonsupersymmetric trilinear scalar interactions, and fi is the 
supersymmetric Higgs(ino) mass, which also enters trilinear scalar vertices. Notice that the 
off-diagonal elements of these squark mass matrices are proportional to the quark mass. 
In the case of the supersymmetric partners of the light quarks mixing between the current 
eigenstates can therefore be neglected. However, mixing between t squarks can be sizable 
and allows one of the mass eigenstates to be much lighter than the top quark. Sbottom mix- 
ing can also be significant if tan/3 3> 1; even in supergravity models with radiative symmetry 
breaking tan/3 can be as large as m t /m b J13J . 



In the presence of squark mixing, the squark-quark-gluino interaction Lagrangian is given 
by 



C m = -iV2g s T a 



Q 



(cos^gi - sin6>qg 2 ) — (sin^gi + cos9 q q 2 ) — - — 



5f a + h.c.,(3) 



where g s is the strong coupling constant and T a are SU(3)c generators. Note that in eq. @ 
we have assumed M- g > 0. 

x We ignore generation mixing between squarks, which in case of the MSSM is only induced radiatively 
by weak interactions. 



2a. Higgs boson decays into light gluinos 

Summing over colors and taking into account the fact that there are two identical particles 
in the final state, the partial decay widths of the CP-even Higgs bosons Hi, H 2 = h and the 
CP-odd boson H 3 = P into a pair of gluinos in the limit m g — > are given by 



T(H k -> gg) = -aM Hk 



")'(&«) 



(4) 



The amplitudes A q k can be written as 



A \,2 = o( s i,2) 9 sin26' g (m 2 q + m 2 g 2 )C (m g ,m g ,mg 2 ) - {ml + m?)C (m g ,m g ,m^) 



~m q sin 26 q 



s ii) Co (m qi , m qi , m q ) - (s 22 2 ) C (m q2 , m q2 , m q ) + 2 (sg) cot 26 q C (m qi , m q2 , m q ) 



A l = -(s 3 ) g sin26 q [(m 2 q - m~ 2 )C (m q ,m q ,m q2 ) - (m 2 q - rri-jCoirriq, m q , m qi ) 



[m qi ,m q2 ,m q ), 



+ (s l2 ) q m q C 
where the scalar function Co(mi,m2,m 3 ) is given in integral form by 

C (mi, m 2 , m 3 ) — — dy dx M 2 ^ x[x — y) + (m 2 — mfjy + (m 2 — m\)x + m 2 — ie 



o jo 



(5) 



"fo 



A complete expression after integration over the Feynman variables is given in the Appendix. 
In the limit Mfj <C m 2 this function simplifies to 



Co{mi,m 2 ,m 3 ) 



m 2 — mi 



m\ log m\ — m 2 log m 2 m 2 log m 2 . — m 2 log m 



2 2 

mf — m 3 



2 2 

m 2 — m 3 



(7) 



and for mi = m 2 



C (m,m,m 3 ) 



mi. 



m\ — m 2 



{m% — m 2 ) 2 



m 2 

lo s — 

mi 



(8) 



The H k qq couplings (sk) 9 are given by || (s 



w 



1 - d = sin 2 6 W ) 



-.2 



m n 



[Sk 



2s\yM W 



fc) 



(9) 



with 



sin a 

sin/3 



cos a 



cos/3 



5 '2 



, r 2 



cos a 
sin/? 

sin a; 

cos/3 



1 



> r 3 



tan/3 
, rn = tan/3. 



(10) 



The Hkqiqj couplings (s^) q are [|TJ 

«')« = MzC ° S(a + ' 3) (/* cos* B, - e A cos 20,) + <f - f^irlA, + r 2 V); 
, ™> s „ M z cos(a + /3) / T r ,„ 9 x 7/1^? m sin 20 . „ . „ . 

( ~22 )9 = _^A ^Z (/£ sin 2 g + ^ CQS 2 ^ + ^1 + 9 {rlAq + r , 2fi) . 

«««=— " -sin2^ g (e q s 2 w --I^ q )-— -*— (rfA, + r^) cos 26 q , (11) 



S«)« = rsi^* 



i 7 ) 9 [sin a — > + cos a , cosa — > — sin a]; (12) 



(~4 1 ) q = (-sf) q = , (§iy = ^-(/i - rlA q ). (13) 

A few remarks are in order: 

(i) As already mentioned, eqs.([5])-(§) are valid only if the gluinos are nearly massless; 
this is an excellent approximation for the case of interest, m- g < 5 GeV <C m t ,mH k ,m q . 
Complete expressions for a finite gluino mass are given in the Appendix. Note that the 
contribution of diagram Fig. lb is ultraviolet finite but the contribution of diagram Fig. la 
is finite only once the contributions of both squarks of a given flavor have been added. 

(ii) The main conventional decay mode of Higgs particles with masses below 130 GeV is 
the decay into bb pairs; in the limit Mjj k 3> m&, the corresponding decay width is given by 

r ( ^5) = ^(^)V!)>. (i4) 

Note that one has to include the QCD corrections to this decay width. The bulk of these 



k ' 



corrections can be absorbed []T5j into running quark masses evaluated at the scale \x = M^ 
For Higgs masses arround 100 GeV, the 6-quark mass mb(mb) = 4.5 GeV will drop to 
the effective value rrib(M Hk ) ~ 3.2 GeV; this results in a decrease of the decay width by 
approximately a factor of two.0 

(Hi) In the limit where either A q and \x or m q are set to zero, there is no mixing between 
left- and right-handed squarks. In this case, the amplitudes A\ are ex rrig and hence very 
small; see also ref. 0]. Therefore only the contribution of the top quark and its SUSY 
partners have to be taken into account; (s)bottom loop contributions are sizable only if 
tan/3 > 1. 

2b. The decay Z — > gg 

The partial decay width of the Z boson into a pair of massless gluinos is given by 



r(^-SS) = ^» W 



4&C 2 W S 2 W V 7T / 






(15) 



2 In contrast to the running of as, light gluinos have very little effect on the running of mb- 



where [here s = Mf ] 

B i = ( v i + a i) a q [ 2 C 2 (s, m q , m q , m q .) - B (s, m q , m q ) - (m 2 + m q 2 )C {s, m q , m q , m qi ) 
+2(viai)v q 2C 2 (s, m q , m q , m qi ) - B (s, m q , m q ) + {m 2 - m q 2 )C (s, m q , m q , m§, 



Bfj = -2a ij {v i a j + Vjai)C 2 (s, m qi ,m qj ,m q ), 



(16) 



with 



an = 2(2/f cos 2 9 q - 2s 2 w e q ) , a 22 = 2(2/f sin 2 6 q - 2s^e g ) , a 12 = -2/f sin 2^, (17) 



2ll L - As 2 e 



w^q i 



21. 



31 

9 ' 



and 



1 1 

V\ = -(cos6 q — sin^g) = a 2 , a± = -(cos 6^ + sin 6^) = — v 2 . 

Alternatively, the result can be written as 



(18) 



(19) 



r(z - ^) 



aMj 



48c%s^ \ ~ i 



asV 



E^ 



with 



5" 



(20) 



= -(a g + i;,cos20,) [2C 2 (s, m q , m q ,m qi ) - B (s,m q ,m q ) + (m q - m qi )C (s,m q ,m q ,m qi ) 
+ -(a q -v q cos29 q ) \2C 2 (s,m q ,m q ,m q2 ) - B (s,m q ,m q ) + (m 2 - m 2 2 )C (s,m q ,m q ,m q2 ) 



- a g m 2 [C (s, m ? , m q , m qi ) + C {s, m q , m q , fn q . 2 )} 

- cos 26 g [auC2(s,m qi ,m qi ,m q ) - a 2 2C 2 (s,m, 2 ,m ( , 2 ,m 9 )] 
+ 2 sin 26 q a 12 C 2 (s, rh qi , m q2 ,m q ). 



(21) 



In terms of the scalar two and three point functions B and Cq given in the Appendix, 
the function C 2 is defined as 



C 2 (s, mi, m 2 , m 3 ) = -m 2 C (s, m 1 ,m 2 , m 3 ) + - + --B (s, mi, m 2 ) - 7r( m i + m l ~ 2m l) 
2B (s, mi, m 2 ) - 5 (0, m 3 , mi) - B (0, m 3 , m 2 ) + (2m 2 - m 2 - m2)C (s, m 1 ,m 2 , m 3 ) 
+^( m i _m 2) 5 o(0, m 3 ,mi) -S (0, m 3 ,m 2 ) + (m 2 - m 2 )C* (s, mi, m 2 ,m 3 ) (22) 

Note that S| and 5?- are ultraviolet divergent [the term proportional to the vectorial 
coupling v q in Bf is finite when one sums over the contributions of the two squarks of a given 
flavor]. It is only when one sums over a complete isodoublet that the amplitudes are finite. 
The contributions from (s)bottom loops therefore always have to be included here. On the 
other hand, if both squarks and quarks of a given generation are degenerate in mass, this 



generation will not contribute to B q \ see also refs. [10 , [XT 



3. Numerical Results 

We are now in a position to present numerical results for branching ratios of Z and Higgs 
bosons into pairs of light gluinos. As stated in the Introduction, we will assume minimal 
(s)particle content, i.e. two Higgs doublets as well as the usual gauge and matter superfields. 
Moreover, unless stated otherwise we will assume that explicitly SUSY breaking contribu- 
tions to the squark mass matrices are identical for all flavors. In terms of the parameters 
introduced in eqs.(|2]), this implies: 

m ~tL = m ~tR = m b R = m V ( 23a ) 

A t = A b = A- rriq. (23b) 

When computing Higgs branching fractions we have to specify one more parameter, besides 
those appearing in eqs.(@). We choose this to be the mass mp of the pseudoscalar Higgs 
boson. We include full 1-loop corrections [0 from the (s)top and (s)bottom sectors to the 



mass and mixing angle of the scalar Higgs bosons, including non-logarithmic terms [[i~7j . 

The values of these parameters are constrained by unsuccessful searches for squarks and 
Higgs bosons. However, many bounds that have been derived under the assumption that 
gluinos are heavy are no longer valid if the gluino mass is just a few GeV. In particular, 
heavy squarks would almost always decay into the corresponding quark plus a gluino. Since 
a light gluino will lose a substantial fraction of its energy in radiation prior to its decay 
[[J, the missing transverse momentum in such events might be too small for them to pass 
cuts designed for conventional squark searches at hadron colliders [|]. On the other hand, 
LEP experiments allow to place bounds on squark masses simply from the measurement of 
the total and hadronic widths of the Z boson, independent of how the squarks decay. We 
therefore require all squark mass eigenstates to be heavier than 45 GeV.Q 

Turning to bounds on the Higgs sector, we note that searches for Z — > Z*h decays (the 
Higgs bremsstrahlung process) are little affected by Higgs decay modes as long as it decays 
hadronically. We have therefore excluded combinations of parameters that violate the m^- 
dependent bound on the ZZh coupling derived by the ALEPH collaboration |18||] On the 



other hand, searches for Z — > hP decays do make use of the assumption that Higgs bosons 
accessible at LEP decay predominantly into pairs of b quarks or r leptons; some sort of 
heavy fermion tagging is necessary in order to suppress the QCD 4-jet background. Thus 
bounds from searches for associate hP production may not apply if these Higgs bosons have 
a significant or even dominant branching ratio into light gluinos. We have therefore not 
included these bounds in our analysis. 

In figs. 2a,b we show results for Br(h — > gg) for rriq = 400 GeV, \x = 200 GeV and several 
combinations of mp and tan/3. The curves have been obtained by varying the A parameter 
in the region A < from its minimal allowed value, defined by mj = 45 GeV, to the point 
where m^ is maximized (At = —\i cot/3) . The results are presented as a function of the light 

3 In practice we only include contributions from third generation (s)quarks; large masses for squarks of 
the first and second generation would not affect our results. This makes it even less likely that the parameter 
space relevant for us is constrained significantly by sparticle searches at pp colliders. 

4 The OPAL collaboration has recently published improved Higgs search limits jl9| . Unfortunately the 
information given in their paper does not allow to extract an upper bound on the ZZh coupling. 



stop mass m- tl , since it is more directly accessible experimentally than the A parameter. 

We find that for mp 3> m% the dependence of the branching ratio on either mp or tan/3 
is quite weak. In this limit the light Higgs scalar h couples to quarks, leptons and massive 
gauge bosons with the same strength as the single physical Higgs boson of the SM P, |6| . If in 



addition ra~ <C ra~, the htiti coupling (s}, l ) q is [1 i\ proportional to the off-diagonal element 



of the stop mass matrix (^a]); this implies that (s^ 1 ) 9 oc (m^ — m~ J /mw oc m~/mw in this 
limit. For fixed mj , the square of this coupling, and hence T(h — > gg), therefore grows like 
the fourth power of rriq. This explains why decays into light gluinos completely dominate 
all other decay modes for small m^ and large mp, as shown in figs. 2. Even for the rather 
modest value of rriq chosen here this new deay mode can suppress the branching ratios into 
bb or t + t~ pairs by more than a factor of 1,000. 

The situation is more complicated for smaller values of mp. If mp is of order rriz or 
less, the scalar mixing angle changes such that T(h —> bb) oc tan 2 /? for large tan/3, while 
the coupling of h to top quarks becomes oc cot/3. Both effects lead to a rapid decrease of 
Br(h — > gg) if mp is not large. Even for tan/3 as small as 2, fig. 2b, for small mp the gluino 
mode is important only if m^ is very small, although it might still remain measurable over 
a wider region of parameter space. Notice also that the branching ratio now depends on A t 
and /i separately, not only on the combination A t + /i cot/? that appears in the stop mass 
matrix (dotted and short dashed curves). 

Fig. 3 shows results for the decay of the pseudoscalar Higgs boson into light gluinos. The 
couplings of P are uniquely determined by tan/3, independent of mp; our result is therefore 
almost independent of mp as long as it lies in the region 1m^ < mp < 2m 4 |] In fig. 3 we 
took mp = 50 GeV, and chose the same values of rriq and /i as in fig. 2. 

We see from fig. 3 that the maximal contribution of the gg final state to the total decay 
width of P is clearly much smaller than in case of the scalar h. The reason is that there 
is no Ptiii coupling M. The Pt\ti coupling (JTJ) can still become quite large if m^ <C rriq; 
however, this coupling only contributes through a diagram containing one heavy stop t 2 in 
the loop, giving a factor l/m~ from the loop integration. As a result, for fixed m^ the 
P — > gg partial width becomes independent of rriq once m? 3> m 2 , as compared to a growth 
oc m| in case of the scalar Higgs boson h. Nevertheless, for small tan/? the gg mode can 
still be quite important even for the pseudoscalar, suppressing the branching ratio into b 
and r pairs by a factor of about 3. For larger tan/? the width for these standard modes 
grows oc tan 2 /3, while the Ptt and Ptrfi couplings decrease; the gg mode therefore becomes 
negligible for tan/3 > 5 or so. 

The same parameters that determine the partial widths for Higgs boson decays into light 
gluinos also determine T(Z —>■ gg); one might therefore hope to place bounds on the former 
by constraining the latter. However, fig. 4 shows that T(Z —>■ gg) depends only very weakly 
on details of L — R squark mixing, described by the parameters A, /x and tan/?. Recall that 
we assume squarks of the first two generations to be degenerate in mass; these generations 
do therefore not contribute to the Z —>■ gg amplitude. Mixing, as well as splitting of mass 



5 We include P — ► Zh decays where kinematically allowed when computing the total decay width of P. 
but this contribution is usually very small. 



eigenstates, can be sizable for stop squarks; we have seen above that the partial widths 
for Higgs decays into light gluinos strongly depend on these details of the stop sector. In 
contrast, the Zt{tj couplings are pure gauge couplings, which are usually modified only 
weakly by squark mixing. Moreover, reducing mi below m t has little effect on the loop 
integrals. It is also important to keep in mind that we get a meaningful result only after 
summing over all four squark mass eigenstates of a given generation; this summation further 
reduces the dependence on the details of the squark mass matrices. Even increasing m~ 
and m~ to 2 ■ m~ (dotted curve) reduces YiZ — > gg) by at most 20%, since the couplings 
of right-handed squarks to the Z are rather weak anyway. 

We see that within our assumptions Br(Z — > gg) < 7 • 10~ 4 . Somewhat larger branching 
ratios are possible Jll| if we allow mass splitting between squarks of the first two generations, 
but this splitting is tightly constrained by bounds on flavor-changing neutral currents, most 
notably K°K° mixing p0[ . The contribution of light gluinos to the total or hadronic decay 
width of the Z boson is thus always well below the present experimental error. f] One will 
therefore have to devise cuts that allow to distinguish gluino pairs from pairs of light quarks. 
One obvious example is to seach for events with missing energy. However, this requires 
knowledge of the gluino fragmentation function which is not very well understood if gluinos 
are light |§. Moreover, there are substantial backgrounds, both instrumental (from the 
mismeasurement of jet energies) and irreducible (from the production and subsequent semi- 
leptonic decay of heavy quarks) . 



A more promising strategy might therefore be |22| to search for g pairs by looking for their 
decay vertices. Farrar recently concluded @ that m g > 3 GeV, mostly due to constraints 
from searches for radiative T decays into gg bound states; the decay vertices of gluinos of 
this mass should be readily detectable, unless the LSP mass is very close to the gluino mass. 
If the LSP is an (almost) massless photino, the gluino lifetime is approximately given by 

io 3 GeV (nig \ , s 

r- q ~ 5 • KT 18 sec — • 24 

m- g \m-gj 

This gives a mean flight path of approximately 0.1 mm (2 cm) for nig = 50 (100) GeV and 
m g = 3 GeV. A detailed Monte Carlo study will be necessary to determine whether detection 
of light gluinos at LEP via these or other |23| , p2| methods is feasible. 

Finally, in figs. 5 a-c we show scatter plots in the planes spanned by any two of the three 
branching ratios discussed above. These plots have been obtained by randomly choosing 
sets of input parameters within the limits 50 GeV < mg < 500 GeV, 50 GeV < \i < 1 TeV, 
25 GeV < m P < 500 GeV, 1.2 < tan/3 < 25 and -4 < A < 4, excluding combinations of 
parameters that lead to too light a squark eigenstate or too large a ZZh coupling. Fig. 5a 
shows a weak positive correlation between the h and P branching ratios: Most points with 
Br(P — > gg) > 10% also have Br(h — > gg) > 1% at least; there are some counter-examples, 
however. Also, the number of points with Br(h —>■ gg) > 10% is clearly much larger than 
that with Br(P — > gg) > 10%; this is not surprising, given the very different dependence of 

6 Light gluinos will also contribute to the hadronic and total Z boson width through loop diagrams; 
however these contributions are generally small fell . 



the two decay amplitudes on tan/? and on rriq, as discussed above. 

Fig. 5b shows that the h and Z branching ratios into light gluinos are anti-correlated. 
The reason is that, as explained above, a large Br(h — > gg) can most easily be obtained if 
m~ 3> m~ is sizable, which suppresses Y[Z — ► gg). Fig. 5c shows no significant correlation 
between the P and Z branching ratios into gluinos. Taken together, figs. 5 a-c demonstrate 
that the measurement of any one of the three branching ratios discussed here would not 
allow one to predict, or even significantly constrain, the other two. 

4. Summary and Conclusions 

In this paper we computed the decay widths of Higgs and Z bosons into pairs of light 
gluinos. Although this decay only occurs at the 1-loop level, it can dominate the total decay 
width of the light scalar Higgs boson by a large factor if the off-diagonal entry of the stop 
mass matrix is large, in which case the lighter stop eigenstate is expected to be (much) 
lighter than the other squarks. The contribution of the gg final state to the total decay 
width of the pseudoscalar Higgs boson can also be quite sizable, although not as large as for 
the sclar Higgs boson. In contrast, the branching ratio for Z — > gg decays is always below 
one permille; a dedicated search will be necessary to explore this possibility experimentally. 
Unfortunately placing bounds on, or a measurement of, Br(Z — > gg) will not teach us much 
about Higgs branching ratios into light gluinos. 

What are the consequences of large Br(h,P —>■ gg)l As discussed in sec. 3, searches 
for Z — > Z*h decays at LEP1, or e + e~ — ► Zh production at LEP2, are probably not much 
affected, since here it is usually not necessary to tag the Higgs decay products in order to 
isolate a detectable signal. Some heavy fermion tagging is necessary in searches for associate 
hP production, however. Notice that at present only the combined limits on Z*h and hP 
production allow one to place lower bounds on rrih and irtp ||TE|, |T5| . We are therefore forced 



to conclude that the present bounds on the masses of the Higgs bosons of the MSSM may 
well be invalid if there are indeed gluinos with mass of only a few GeV. 

On the other hand, in the long run the existence of light gluinos need not hamper Higgs 
searches at e + e~ colliders. While gluino pairs do look different from b and especially r pairs, 
they also ought to differ sufficiently from standard light quark or gluon jets to allow the 
suppression of QCD backgrounds. The impact of light gluinos on Higgs searches at hadron 
colliders might be less benign, however. The huge background from pure QCD processes 
means that gluino pairs cannot be used as a Higgs signal here. We have seen that in the 
presence of light gluinos the branching ratio of scalar Higgs bosons into b or r pairs might 
be reduced by a factor of 1,000 or more. The same suppression factor applies for decays 
into (virtual) gauge bosons. The h — > ZZ* — > 4 leptons signal would then be undetectable; 



however, within the MSSM it is at best marginal anyway ||24j| . 

Potentially most unfortunate for the prospects of MSSM Higgs searches at hadron col- 
liders would be a reduction of the h — > 77 signal, which is usually regarded to be the most 



promising way to search for h [^4] . However, the same ht\t\ coupling that can give rise to a 



large partial width for the gg final state also contributes to the h — > 77 and h — > gg partial 
widths, via t\ loops. A proper investigation of the possible impact of the existence of light 
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gluinos on Higgs boson searches at hadron colliders therefore necessitates an analysis of all 
SUSY loop contributions to Higgs production and decay, which is beyond the scope of this 
paper. We nevertheless hope that our results will give additional urgency to the effort to 
either detect or definetely exclude the existence of gluinos with mass of a few GeV. 
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Appendix: Complete result 

Taking into account the mass of the gluino, the partial decay width of the CP-even and 
CP-odd Higgs bosons into gluino pairs is given by 



T(H k -► gg) = ^M Hk 



1 .. /«s\ 2 Rl , 



IT J 



1% 



»J=1 



.1=1 



(A.1) 



where f3 g — 1 — AM~/M]j is the velocity of the final gluinos with p — 1 for fc = 1, 2 and 
p = 3 for fc = 3. The amplitudes A\ come from the diagrams Fig. la and A % £ come from the 
diagrams Fig. lb; they are given by 



^1,2 = (si t2 ) q (Vi - a { ) [B (s, m q , m q ) + (m g +m q + m,- )C (s, m q , m q , m qi ) + Am g 

C+(s, m q , m q , m q .)} + 2m q m g (v 2 + a 2 ) [C (s, m„ ra 9 , m$J + 2C+(s, m 9 , m„ m^)] 



^3 = ( s 3) q ( v i ~ a i ) ~ B o(s, m q , m q ) + (m g + m - m$. )C (s, m„ m„ to,.) 



(A.2) 



+2m q m- g (v 2 + a 2 )C (s, m q , m q , m qi ) 
A % 1 2 = (si, 2 ) m q (viVj - OiOj)C (s, m qt ,m qj ,m q ) - 2m g (viVj + aiaj)C + (s, m qt ,m qj ,m q 
m q (via,j - Vjai)C (s, m qv m qj ,m q ) - 2m g (v i a j + aj^)C_(s, m^, m qj ,m q ) 



A'i 



with the couplings (sfc) 9 and (s^) 9 given in eqs. (6-10) and 

1 1 

V\ = -(cos9 q — smO q ) = a 2 , <X\ = -(cos 6^ + sin# g ) 

The functions C + (s,mi,m 2 ,m 3 ) is defined as 

1 



-v 2 



(A.3) 



C + (s,mi,m 2 ,m 3 ) 



2^1 
1 



2B (s,m 1: m 2 ) - B (m g 2 , m 3 , mi) - B (m g 2 , m 3 , m 2 ) 
+ {2m g 2 + 2m 2 — m 2 — m 2 )Co(s, mi,m 2 , m 3 ) 



(A.4) 



C_(s,mi,m 2 ,m 3 ) = - 1 [So(mg 2 ,m 3 ,m 2 ) - B (m g 2 ,ra 3 ,mi) + (m\ - m 2 ) 2 C (s,m 1 ,m 2 ,m 3 ) 



with the scalar two and three point functions, B and Cq defined as 

(2n/j,) n - 4 r d n k 



B Q (s,m 1 ,m 2 ) 



C (s,m 1 ,m 2 ,m 3 ) 



27T Z 



(27T/X) 



n-4 



11X Z 



(k 2 — m 2 + ie)[(k — q) 2 — m 2 + ie] 

d n k 



(A.5) 



[(fc — pi) 2 — m 2 + ie] [(k — p 2 ) 2 — m 2 + ie] (k 2 — m 2 + ie) 



where n is the space-time dimension and \x the renormalisation scale. After integration over 
the internal momentum k, the function Bq is given by [je is Euler's constant] 



„ , ,1 , 4iru 2 „ m? — m 2 . . m 2 . x+ — x_ 

S (s,mi,m 2 ) = --7s + + log — — + 2 + ^- -log-| + -^- 

e mi7n 2 2s m{ As 



1 x " 
log — 

X_|_ 
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with 



x± — s — ml — m\ ± J s 2 — 2s(m 2 + m 2 .) + (m\ — m 2 .) 2 



(A.6) 



1 f y r 2 2 i _1 

ay / dx ay + bx + cxy + dy + ex + / 

o jo L ' ' J 



(A.7) 



The three point scalar function Co for p\ = p\ = m~ and s = (p\ + P2) is given in integral 
form by 

C ( y s,mg,m 1 ,m 2 ,m 3 ) = 

where 

a = rrig 2 , b — s , c = — s , <i = m2 — mj — m^ 2 , e = m 2 — m 2 , , / = m 2 — ie. (A. 8) 
Co can be expressed in terms of a sum of Spence functions Li2(x) = — J dt\og(l — xt)/t: 

Xi - 1 



C (s,m g ,m 1 ,m 2 ,m 3 ) 
where we have defined 






Li? 



,/.. ; 



v <Ei Vij , 



Li 2 



v •t'l i/jj / _ 



,(A.9) 



Xi 



x 2 



x 3 



2d + e(l-(3 q ) 1 

2rf + e(l-/? g ) 
s/5,(l + /3,) 
2rf + e(l-fl,) 
sp q (l-p q ) 



P, 



1) 1 



V3± 



— c — e ± w(c + e) 2 - 


- Ab(a + d + f) 


2b 


-d-e±^J(d + e) 2 


-4f(a + b + c) 



2{a + b + c) 
-d ± y^ 2 - 4a/ 
2a ' 



(A.10) 
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Figure Captions 



Fig.l The decays Z,h,P —>■ gg proceed via 1-loop diagrams with two quark propagators and 
one squark propagator (la), or two (possibly different) squark propagators and one 
quark propagator (lb). 



Fig.2 The branching ratio for the decay of the light scalar Higgs boson h of the MSSM into 
a pair of light gluinos is shown as a function of the mass of the lighter stop eigenstate 
ti, for the ratio of vacuum expectation values tan/? = 25 (2a) and 2 (2b), respectively. 
We have fixed rriq = 400 GeV, and ii = 200 GeV for all cases except the dotted curve 
in fig. 2b, which is for /i = 500 GeV. The curves have been obtained by varying the A 
parameter in the region A < 0, as explained in the text. 



Fig.3 The branching ratio for the decay of the pseudoscalar Higgs boson P of the MSSM 
into a pair of light gluinos is shown as a function of mi , for rriq = 400 GeV, /x = 200 
GeV, mp = 50 GeV and three different values of tan/3 as indicated. The curves for 
small tan/? terminate well below the possible maximum of m~ tl since here substantial 
1-loop corrections to the scalar Higgs sector are necessary to evade the ALEPH bound 
|18| on the ZZh coupling; this leads to a lower bound on \A t \ in these cases. 



Fig.4 The branching ratio for the decay of the Z boson into a pair of light gluinos is shown 
as a function of the mass rriq L of SU(2) doublet squarks, for /i = 200 GeV and various 
combinations of A and tan/3. The solid and dashed curves have been derived using 
our usual assumption of equal explicit SUSY breaking masses for SU(2) doublet and 
singlet squarks, while the dotted curve is for m~ fl = 2m ~ L . 



Fig.5 Scatter plots in the planes spanned by the branching ratios of h and P (5a), h and Z 
(5b) and P and Z (5c) into pairs of light gluinos. These plots have been obtained by 
randomly choosing combinations of input parameters within the boundaries 50 GeV 
< m~ q < 500 GeV, 50 GeV < // < 1 TeV, 25 GeV < m P < 500 GeV, 1.2 < tan/3 < 
25 and —4 < A < 4. Combinations of parameters that violate the bounds on squark 
masses and the ZZh coupling discussed in the text have been discarded. 
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